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ABSTRACT 


The intermittent behavior of solar wind turbulent fluctuations has often been 


investigated through the modeling of their probability dis 


(PDFs). Among others, the Castaing model fjCastaing et ah 


ribution functions 


1990[) has success¬ 


fully been used in the past. In this paper, the energy dissipation field of solar 
wind turbulence has been studied for fast, slow and polar wind samples recorded 
by Helios 2 and Ulysses spacecraft. The statistical description of the dissipation 
rate has then be used to remove intermittency through conditioning of the PDFs. 
Based on such observation, a self-consistent, parameter-free Castaing model is 
presented. The self-consistent model is tested against experimental PDFs, show¬ 
ing good agreement and supporting the picture of a multifractal energy cascade 
at the origin of solar wind intermittency. 


Subject headings: solar wind, turbulence, intermittency 
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Introduction 


The interplanetary space is permeated by the solar wind, a rarehed, magnetized 
plasma continnously expanding from the solar corona. The solar wind blows radially 
away from the snn, and extends np to abont lOOAU, at snpersonic and snperalfvenic 


speed. Measnrements collected by spacecraft instrnments dnring 


evidenced that low freqnency flnctnations have power law spectra (jColeman 


ast decades have 


19681 1. This 


snpports th e stndy of flnctnations in the framework of magnetohydrodynamic (MHD) 


tnrbnlence flBruno and Carbone 


2013 


Biskamp 


1997 1. The MHD natnre of the tnrbnlent 


flnctnations has recently been c onhrmed by more detaUed analysis of the linear scaling o f 


the mixed third ore 

er moment 

Politano anc 

PouQuet 

1998; 

MacBride et al. 

Sorriso-Valvo et al. 

2007, 

201 ol 

Marino et al. 

2008, 

2011, 

2012; 

Carbone et al. 


2005 


2008 


200911 . One 


of the most interest ing properti es of solar wind tnrbnlence is the intermittent character 


of the flnctnations fiFrisch 


199511 of helds snch as velocity, magnetic held, or the Elsasser 


helds. Intermittency is related to the non-homogeneons generation of energetic strnctnres 
in the how, due t o the nonlinear transfer of energy across the scales a s observed in 


geophysic al hows llPonanet et al 


plasmas llBrnno and Carbone 


cross-helicity flSmith et al. 


2013 


Marino et al. 


2013a . 


2014 . 


20151 ) and heliospheric 


20131 ). whose efficiency can be corre 


200911 or self-generated kinetic helicity flMarino et al.ll2013b|l . 


ated with levels of 


Being ubiquitous in turbulence, intermittency plays a relevant role in the statistical 
description of the held huctuations. The main manifestation of intermittency in fully 
developed turbulence is the scale dependent variation of the statistical properties of the 
held increments, customarily dehned as = 'ip{x + £) — 'ip{x) for a unidimensional 

generic held V’(^) at the scale £. In particular, many studies have focused on scaling 
properties of the probability distribution functions (PDFs) of the held huctuations, P{A'ipi), 
and of the structure functions, dehned as the moments of the distribution function of 
the held huctuations Sp{A'ijji) = {Atf^i^xy) (here (■) indicates an ensemble average). The 
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intermittent (i.e. non-homogeneous) concentration of turbulent energy on small scale 
structures, as the energy is transferred through the scales, results in the enhancement of 
the PDFs tails, indicating that large amplitude fluctuations are more and more probable 


at smaller and smaller scale. Correspondingly, the structure functi ons sea ' 


deviate from the linear prescription valid for scale invariant PDF (iFrisch 


ing e xponents 


19951 ). Such 


deviation has been studied since early wor ks on turbulence. Metho ds based on wave let 


transform flFarge 


1992 : 


Onorato et ah 


20001) or threshold techniques flGreco et al.ll2009f) have 


been developed for the identiheation, description and characterization of the intermittent 


structures. On the ot 


PDFs fICastaing et al 


rer hand, models for intermittency try to reproduce the s. 


1990 


Benzi et al. 


the structure functions anomaly flMeneveau and Sreenivasan 


1991 


Frisch and Sornette 


19871 : 


19971: 


rape o f the 


Frisch 


1995h or 


She and Leveque 


19941) . 


allowing the quantitative evaluation of intermittency. 


In solar wind plasma, interm i ttency has been extensively characetrized in th e 


recent years flTu and Marsch 


1995 


Marsch and Tu 


1997- 


Sorriso-Valvo et al 


1999 


2001 


Bruno et al. 


20091: 


2007a: 


Kivani et ah 


Hnat et al. 


20091: 


Ragot 


2003 


Lenbner and Voros 


2005a; 


Macek and Wawrzaszek 


20131) . As result of intermittent processes, small scale 


structures, like for exa mple thin current shee 


identih ed in solar wind fjVeltri and Mangenev 


:s or tangential discontinuities, have 

3een 

1999; 

Greco et al. 

2009; 

Bruno et al. 

O 

O 

_ 


2007bl) . Their dependence on parameters such as solar 


activity leveh heliocentric distance. 


heliolatitude and wind 

speed has been a 

so pointed out ( 

Pagel and Balogh 2002c 

Bruno et al. 

0^ 

o 

o 

CM 

Nicol et al. 

20Qa; 

Ghapman et al. 

200§; 

Yordanova et al. 

2009 

)• 


One of the mode 


s for the description of the increments PDFs was introduced by Gas- 


taing 

Gastaing et al. 

1990 

), anc 

successfully applied in several contexts ( 

Sorriso-Valvo et al. 

1999, 

2000, 

2001 

Padhve et al. 

200 


Page] 

and Balogh 

2002 

: Stepanova et al. 

2003 


Hnat et al. 

2 OO 3 I; 

Lenbner and Voros 

200,5b, f 

n Garbone et al. / 

>004 

). Within the multifractal 
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framework flFrisch 


1995f) . for each scale, the energy transfer rate has non-homogeneons 


scaling properties, as for example the fractal dimension (which is directly related to the 
cascade efficiency) in different regions of space. The tnrbulent helds can be thns interpreted 
as a superposition of subsets, each characterized by a given fractal dimension, and with a 
typical energy transfer rate. Each region can then be reasonably assumed to have the same 
distribution Pq of the held huctuations, with variable width a (depending on the cascade 
efficiency, and related with the local fractal dimension) and weight L{a) (depending on the 
fraction of space characterized by the same statistics). The Castaing model PDF consists 
thus of the continuous superposition of such distributions, each contributing to the statistics 


with its appr opriate weight. The 


the widths a (ICarbone et ah 


atter is introduced through the distribution function of 


20041 ). This leads, for each time scale i, to the convolution: 


F(A^/>£) = / Lxia^)Po{A^jJe,ae) d \nai. 


( 1 ) 


Based on empirical large scale PDF shape, a Gaussian parent distribution is normally 
used. It is known from turbulence studies that PDFs of fluctuations have to be skewed. 



to account for this, a skewness parameter a* must also be included, so that 


Po{A'il:i,ai,as) = 


V2 


: exp 


ira 


M 

2a2 


1 + 


Atl^i/a 


(1 + Atlj^a^yP 


The PDF of variances Lx{(T() needs theoretical prescription. A log-norma 


( 2 ) 

ansatz has been 


ofteii used, as conjectured in the framework of the multifractal cascade fICastaing et al 
199011 


2A| 






exp 


(3) 
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Such choice has been justihed by assuming that the nonlinear energy transfer is the result 
of a multifractal fragmentation process, giving rise to a multiplicative hierarchy of energetic 
structures. By assuming random distribution of the multipliers (namely, of the local 
efficiency of the cascade), the central limit theorem suggests a log-normal distribution of the 
local energy transfer rate. Then, from dimensional considerations, the fluctuations variance 
can be expected to share the same stati stical properties as th e energy transfer rate, therefore 


giving the log-normal distribution (|3]) flCastaing et al 


199011 . In equation 


), for A» = 0, 


the log-normal PDF is a 5-function, so that the convolution ([T]) gives one Gaussian of width 
ao/ the most probable value of As increases, the convolution includes more and 
more values of and the PDF tails are enhanced. Therefore, the scaling of the parameter 
A^ controls the shape of the PDF tails, and describes the deviatio n from the parent 


distribution, charac 


Sorriso-Valvo et al 


observed for A^ fICastaing et al 


;erizi ng the intermittency in the inertial range (jCastaing et al. 


1990 


199911. In fully d e veloped turbulence, a po wer-law scaling is usually 


19901: 


Sorriso-Valvo et al. 


19991) . Finally, a relationship can 


be es 


ablished between th e scaling exponent of A^ and the multifractal properties of the 


flow (jCastaing et al. 


19901 1. 


As briefly described above, the Castaing model is based on hypotheses on the physical 
processes governing the turbulent cascade. The main assumption is the choice of the weights 
distribution, which needs appropriate theoretical modeling. In this paper, after verifying 
that the multifractal cascade framework applies to solar wind turbulence, we show that it is 


possible to describe solar wind intermittency wit 


distribution, but rather using empirical weights (INaert et al. 


rout any hypothe sis on the shape of such 


199811 . In Section [2] we briefly 


introduce the data used for the analysis; in Section [3] we estimate the empirical distribution 
function of the local energy dissipation rate; Section |4] shows the conditioned analysis 
performed on the data in order to extract the empirical weights for the Castaing model, the 
self-consistent Castaing probability distribution functions are built and compared with the 
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experimental PDFs. 


2. The data: Helios 2 and Ulysses 


This work is based on the analysis of three different samples of in situ measurements 
of velocity v, mass density p, estimated using proton and a particle, and magnetic held 
B. The Elsasser variables = v ± B/^/dvrp have also been evaluated from the time 
series. Two samples were taken in the ecliptic wind by the Helios 2 spacecraft during 
the hrst 4 months of 1976, when the spacecraft orbited from lAU on day 17, to 0.29AU 
on day 108. Data resolution is 81 seconds, and eleven fast or slow wind streams, each 
about 49 hours long, 
a better stationarity 


lave been extractec 

to avoid stream-stream interfaces and to ensure 

Sorriso-Valvo et ah 

1999 

). As pointed out in the literature, fast 


and slow wind turbulence shou 


conditions f|Rruno and Carbone 


d be studied separately, because of the different plasma 


20131) . Therefore, for our statistical analysis we built two 


distinct samples by putting together six streams for the fast wind (which we name the fast 
sample, totalizing = 12288 data points), and hve streams for the slow wind (hereafter 
the slow sample, consisting of = 10240 data points). The third sample was recorded in 
the solar wind out of the ecliptic, by instruments on-board Ulysses spacecraft, during the 
first eight months of 1996. The spacecraft spanned distances in the range from 3.1AU to 
4.7AU, and latitudes from about 53° to 19°. Sampling resolution is 8 minutes. We refer 
to this dataset as polar sample. All samples were taken at solar minima, when the solar 
wind is more steady, and free from disturbances of solar origin. We remind that, in order 
to study spacecraft time series, all scales i are customarily transformed in the time lags 
T = f'/|(v)| through the bulk flow speed averaged over the entire data set. This is allowed 


by the 


Taylor hypothesis, tha 


range fjPerri and Ralogh 


2010 


is generally valid for solar wind fluctuations in the inertial 


Rruno and Carbone 


20131 ). and which we have tested in 






















our data. In Figured], samples of the PDFs of standardized magneti c and velocity field 


fluctuation s at different scales are shown for fast arid pola r samples fjSorriso-Valvo et ah 


1999, 

2001 ; 

Carbone et al. 

2004; 


Sorriso-Valvo et al 


toward small scales appears evident for all cases (ISorriso-Valvo et ah 


20041) ■ The typical ta il enhancement 


are o 


as erved for the yelocity and for the Elsasser helds (see e.g. 


lllOOOh : 


19991. Similar results 


Sorriso-Valvo et al 


Carbone et al. 


(1200411 1. PDFs of magnetic held, velocity and Elsasser helds 


huctuations have been successful 


Equation (dj) fISorriso-Valvo et al. 


y rep roduced using the Castaing model given in 


19991) . In particular, the heavy tails are extremely well 


captured by the model, showing that the ehect of intermittency has been well described. 
The htting parameter shows a power-law scaling extended over about one decade. As 


mentioned above, the scaling exponents are related to the fracta 


intermittent structures generated at the bottom of the cascade fjCastaing et al. 


dimension of 


le most 


199011 . 


therefore providing important physical information on the turbulent how properties. 


3. Statistical properties of the energy dissipation rate 

The hrst step for the self-consistent characterization of the PDFs is to map the 
properties of the plasma energy transfer rate using the experimental time series. In solar 
wind turbulence, understanding the mechanism responsible for the dissipation (and/or 
dispersion) of the energy at the bottom of the nonlinear cascade is still an open issue. 
Therefore, the actual expression of the dissipative (and/or dispersive) terms within MHD 
equations is not known. It should also be pointed out that, unlike ordinary turbulent 
hows, the solar wind plasma is only weakly collisional, so that molecular viscosity and 
resistivity cannot be dehned in a simple way, nor estimated directly from the measurements. 
These considerations show that it is not possible, to date, to measure the local energy 
dissipation rate in solar wind turbulence. However, it is possible to dehne proxies of the 
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T = 81 sec T = 1 hrs x = 11 hrs 



X = 8 min x = 1 hrs x = 1 hrs 



-4-3-2-1 0 1 2 3 4 -4-3-2-1 0 1 2 3 4 -4-3-2-1 0 1 2 3 4 


5bj. 6bj. 5bj. 


Fig. 1.— The PDFs of standardized field fluctuations at three different scales (increasing 
from left to right), for fast solar wind velocity fluctuations from Helios 2 (top panels), and 
for polar Ulysses magnetic fluctuations (bottom panels). 
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energy dissipation rate, that can be reasonably used to represent the statistical properties 
of the held. I n this paper, we use a dehni t ion based on the third or der moment scaling 


law for MHD fjPolitano and Pouquet 


1998 


Sorriso-Valvo et al. 


20071) . often referred to as 


Politano-Pouquet law (PP). The PP law establishes, under given hypotheses (stationarity, 
homogeneity, isotropy, incompressibility), the linear scaling of the mixed third order 
moment of the Elsasser helds. 


Y^{t) = (|Az^,(t)p Az^(t)) = -- (e±)r(n) . 


( 4 ) 


In the right hand side of Equation (jlj), (e^) is the mean energy transfer rate, estimated over 
the whole domain. By analogy, we dehne the “local” pseudo-energy transfer rate as: 


= 


t{v) 


( 5 ) 


so that the local energy transfer rate at the scale r reads e,-(f) = + e^{t))/2. 


At a given scale, eac 


1 helc 


er{t) ( Marsch and Tn 


increment can thus be associated with the local value of 


19971 ). Since we are interested, in particular, in the small scale 


intermittent effects, from now on we will only use the resolution scale values of e{t) (namely 
81 seconds for Helios 2 data and 8 minutes for Ulysses data), which we will refer to as 
energy dissipation rate. Figure [2] shows the variable e{t), computed for fast and slow 
Helios 2 streams and for the Ulysses polar sample. Differences between the three samples 
are evident, in particular for the Ulysses dataset, probably because of the lower data 
resolution. For all cases, the field is highly irregular and inhomogeneous, with spikes of 
large dissipation alternated with quiet periods. The probability distribution functions of 
e(t) are shown in Figure [3] for fast and slow wind, as obtained from Helios 2 data, and 
for polar Ulysses data. Because of the inhomogeneity, PDFs have been computed using 
bins of variable width, by imposing a fixed number of data points in each bin {Np = 1500 
for Helios datasets, Np = 3000 for Ulysses data). The error bars are estimated as the 
counting (Poisson) error on Np, and the three PDFs have been vertically shifted for clarity. 
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Local dissipation e(t) [10^ x J Kg / sec] 



0 2000 4000 6000 8000 10000 


Time index 


Fig. 2.— The local energy dissipation rate (see text) for the three datasets: Ulysses fast 
polar wind (top), Helios 2 fast streams (center), and Helios 2 slow streams (bottom). 
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Log-Normal fits of the distributions are overplot ted in light grey, sho wing good agreement 


with the framework of a multiplicative cascade fjCastaing et ah 


1990 ). Alternatively, PDFs 


are even better reproduced (with two to ten times smaller by a stretched exponential 
£t P(e) ~ exp —(e/eo)'^, where eo is a characteristic value of the energy dissipation rate, 
and c is the parameter controlling the shape of the tails of the PDF. In particular, c = 2 
indicates Gaussian, c = 1 exponential, and c < 1 indicates heavy tailed, almost power-law 
distribution. In the present case, the shape parameters are indicated in Figure |3l Both 
fast and slow wind measured by Helios 2 show a strong deviation from Gaussian (c ~ 0.2), 
with presence of energy dissipation bursts that result in the heavy tails of the PDFs. 

On the contrary, Ulysses shows “thinner” tails, i.e. populated with less extreme events. 
This is to be expected, because of the different time resolution of the data. Indeed, while 
Ulysses sampling time sits in the middle of the inertial range, for Helios 2 it is closer to its 
bottom. It should be recalled that stretched exponential PDFs were predicted for a variable 
generated as the result of a multiplicative pr ocess controlled by a few e xtreme events, in the 


framework of the extreme deviation theory (iFrisch and Sornette 


19971 ). This supports once 


more the presence of a intermittent, multiplicative energy cascade in solar wind turbulence. 


4. The conditioned analysis and a self-consistent Castaing model 

In order to verify that the Gastaing model can properly be applied to solar wind 
turbulence, a conditioned analysis of the data was performed. The range of values of the 
energy dissipation rate was divided into Nbin bins of variable width (with Nbin = 8, 7,10 
for fast, slow and polar samples, respectively), in order to separate different regions of the 
time series. The PDFs of held huctuations, at the resolution time lag, were thus estimated 
separately for each bin, i.e. conditioned to the values of the energy dissipation rate. Top-left 
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Fig. 3.— The probability distribution function of the energy dissipation rate e estimated 
from the Helios 2 data for fast wind (blue circles) and slow wind (red squares), and from the 
Ulysses polar data (green triangles). Fits with log-normal (grey lines) and with a stretched 
exponential function (see legend) are superposed, and the corresponding shape exponents c 
are indicated. The three distributions have been vertically shifted for clarity. 
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panel of Figure H] shows, superimposed, all the ten conditioned PDFs P{A'ijj\e) of the 
resolution scale radial velocity held huctuations for the polar wind dataset. Conditioned 
PDFs no longer have the heavy tails observed at small time lags (see Figured]), and present 
roughly Gaussian shape. The top-right panel of the same Figure shows examples of the 
Gaussian hts of the conditioned distributions. For clarity, only four out of ten conditioning 
values of e are plotted. This observation conhrms the multifractal scenario, in which regions 
with the same energy dissipation rate are characterized by self-similar (Gaussian in this 
case) PDFs of the held huctuations, even at small scale. Thus, the conditioned data lose 
intermittency, which indeed arises from the local huctuations of the energy dissipation rate. 
Top panels of Figure ID also qualitatively highlights the result of superimposing several 
Gaussian distributions of variable width and amplitude, obtaining the heavy tails observed 
in the usual PDF. The use of the Gastaing distribution is therefore appropriate to describe 
solar wind intermittent turbulence. For each conditioned PDF, corresponding to a given set 
of constant dissipation rate, it is possible to obtain the probability density function of the 
standard deviations L{a) = Np/Ntot^e, where Np is the (hxed) number of points in each 
dissipation bin. This is done for each dataset by exploiting the correspondence between a 
given bin of energy dissipation rate e and the standard deviation of the corresponding held 
huctuations. In other words, for each held V’ and for each value of e, a value of cr(A'^|e) is 
obtained through the Gaussian ht of the conditioned PDFs P{A'ijj\e). Then, the probability 
L{a) is the fraction of data associated with that interval of e values. This allows to establish 
an explicit relationship between the statistical properties of the energy dissipation rate e and 
of the standard deviation a, postulated in the Gastaing model. Gentral panels of Figure 0] 
shows the probability distribution functions of the standard deviation of the helds L{a) 
associated to the distributions of e, along with log-normal and stretched exponential hts. 
In this case, the log-normal ht is not always able to capture the form of the distribution, 
while the stretched exponential hts are very good. For the latter, the smaller exponents c 
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Fig. 4.— Top-left panel: conditioned PDFs P(A'0|e) for the radial component of the velocity 
in the polar Ulysses sample, for ten different valnes of e. Top-right panel: the Gaussian £t of 
F(A'^|e) for the same case, only shown for four values of e. Central panels: PDFs of a{A'ip\e) 
for the three dataset, for the radial component of velocity (center-left panel) and magnetic 
field (center-right panel). Fits with a stretched exponential law are superposed (solid lines), 
and the htting parameters c are indicated. Bottom panels: relationship between the energy 
dissipation rate e and the variance a{A'ijj\e) for velocity (bottom-left panel) and magnetic 
held (bottom-right panel) huctuations. 
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observed for the velocity fluctuations in all dataset, and suggesting a lower intermitten cy, 


agree with the higher level of intermittency of the magnetic held fjSorriso-Valvo et al 


19991). It is also interesting to check directly from the data the validity of the assumption 


originally used in the Castaing model, namely the phenomenological relationship between 
the energy dissipation rate and the variance of the huctuations. This can be simply done 
by plotting e versus a{A'ilj\e), and verifying the power-law relation e ~ cr" expected from a 
simple di mensional argument on the third orer sc aling relation given by Eq. (jl]), resulting 


in a = 3 fjCastaing et ah 


1990 


Naert et al. 


19981) . As can be seen in the bottom panels 


of Figure IH the power-law relationship is very well verihed, with exponents compatible 
with the expected value a = 3. The error bars are estimated as half the size of the energy 
dissipation rate bin. Similar results apply to the other components of the helds and to 


the Elsasser helds (not shown). No 


:e that a similar analysis performed on ordinary huid 


turbulence provided instead a = 5 flNaert et al. 


1998|). 


With all this information at hands, the self-consistent Castaing distribution ([T]) can be 
hnally discretized as follows: 


^bin 

P^A'lpr) = Anorm ^ T((T)Po(A?/’r, CT, a^) A,. (6) 

i=l 

Equation ([6]) only includes two free parameters: the normalization factor Anorm, and the 
skewness factor a^, while the information on the PDF shape is self-consistently included 
through the distribution L(a). The PDFs of the helds huctuations can now be compared 
with the self-consistent Castaing model, by htting the data with Equation ([6]) in order to 
adjust only for the two free parameters. These are however irrelevant for the tails curvature, 
i.e. for intermittency, so that there is no need to discuss them here. Figure |5] shows 
the PDFs of huctuations, as in previous Figure [U (markers), together with the ht with 
Equation ([6]) (full line). The self-consistent model reproduces the shape of the PDF tails 
very satisfactorily in all cases. In order to test the goodness of the ht, a synthetic Lji{a) 
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has been built by averaging over 10^ realizations of random distributions L{a), generated 
for a in the same interval as each solar wind dataset, [(Tmin, ■■■,crmax]- The fit of the data 
with the self-consistent Castaing model, using the synthetic Lji{a), is shown in Figure E] as 
dotted line. It is evident that a superposition of Gaussians with a random choice of the 
weigth distribution does not allow the correct reproduction of the data. Quantitatively, the 
random Gaussian superposition fits have values larger by a factor of ten with respect to 
the self-consistent Gastaing £t. 


5. Summary 


In the attempt to understand the mechanism responsible for the generation of the 
highly intermittent turbulent fluctuations of solar wind helds, we have explored here the 
existence of a non-homogeneous energy cascade, a base hypothesis of several models. The 
intermittent solar wind held huctuations, deeply investigated in the last decade, have 
often been described through the Gastaing PDF model. In this work, we have evaluated 
empirically the weights distribution L{a) of the Gastaing model (Eq. ([1])), rather than 
assuming an analytical prescription (usually a log-normal distribution). For the sake of 
generality, we have selected three different datasets, consisting of one fast and one slow 
ecliptic wind samples (Helios 2 data), and one fast polar wind sample (Ulysses data). 
Since it is not possible, to date, to measure the energy dissipation rate from the data, we 
have used a proxy derived from the exact statistical scaling law for the mixed third order 
moment of the magnetohydrodynamic huctuations. Although such law is expected to apply 
only for homogeneous, isotropic, incompressible MHD turbulence, its validity in the solar 


wind time series has recently been established by several authors (e.g. 


Sorriso-Valvo et ah 


fll999l) . Therefore, the use of the proxy introduced here is appropriate. After estimating the 
local value of the energy dissipation rate, we have evaluated its statistical properties. The 
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HELIOS 2 - fast ecliptic solar wind 



5z-j. 


HELIOS 2 - slow ecliptic solar wind 



5z+j. 


ULYSSES - fast polar solar wind 

10 


(O 0.1 


0.01 


0.001 

-4 -3 -2 -1 0 1 2 3 4 

8bj. 


Self-consistent Castaing 
Random Castaing 



Fig. 5.— The PDFs of radial Elsasser or magnetic fields fluctuations at the resolution scale, 
computed for the three samples, together with the self-consistent Castaing model (full line). 
The dotted line indicates the model PDF obtained for random distribution of the variances 
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probability distribution functions are consistent with a log-normal function, predicted by 
the typical theoretical models used for the description of the turbulent helds as the result of 
a non-homogeneous, random multiplicative cascade. Alternatively, a stretched exponential 
£t was also used for all datasets. Stretched exponentia 


framework of the Extreme Deviations Theory (EDT) ( 


distributions mate 

1 the theoretical 

Frisch and Sornette 

1997 

), where 


the statistical properties of the multiplicative cascade process are controlled by the most 
extreme intermittent events (the heavy tails of the PDFs). The stretching parameter, 
obtained from the £t, can be used to describe quantitatively the degree of inhomogeneity, 
being related to presence of correlations or clustering in the dissipation held. The results 
obtained here (see Figu re [3D confirm the stronge r intermittency of slow win d with respect 


to fast and polar wind (jSorriso-Valvo et ah 


1999 


Rruno aud Carbone 


2ni3h . 


Starting from this observation, we have then studied the conditioned PDFs of the fields 
huctuation, the conditioning parameter being the local energy dissipation rate. As expected 
within the multifractal picture (but never observed in the solar wind so far), conditioning 
results in loss of intermittency, so that self-similarity is restored when the inhomogeneity 
of the energy dissipation rate is removed. This result shows that the framework of the 
multifractal energy cascade applies to solar wind turbulence. 


Upon observation of their Gaussian shape, the conditioned PDFs have hnally been 
used to build the Castaing distribution self-consistently, through the empirical distribution 
of their standard deviations, L{a). The model PDF drawn following this procedure hts very 
well the shape of the experimental PDFs, and in particular can capture the curvature of the 
tails. The self-consistent Castaing model can therefore reproduce the intermittency of the 
helds. We remark that, the selection process for the conditioned PDFs being based on the 
dissipation rate, the relevant physical information necessary to describe the intermittent 
PDFs is the non-homogeneous distribution of the local energy dissipation rate. The role of 
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the shape parameter A^, used to describe quantitatively the degree of intermittency, is now 
played by the stretching parameter c of the distributions of the energy dissipation rate, not 
directly involved in the self-consistent Castaing model. 

In conclusion, the stretched exponential distribution of the local energy dissipation 
rate, the Gaussian shape of the conditioned PDFs, and the successful application of the 
self-consistent Castaing model, all strongly support the picture of the multifractal energy 
cascade being at the origin of the intermittency of solar wind turbulence. The validity of 
such scenario holds for both fast, slow and polar samples, confirming that the turbulent 
cascade is always active in the solar wind. 

The research leading to these results has received funding from the European 
Communitys Seventh Framework Programme (FP7/2007-2013) under grant agreement n. 
269297/ “TURBOPLASMAS”. 
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Helios2-R=0.9 AU 































































































